Electronic structure via potential functional approximations 
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The universal functional of Hohenberg-Kohn is given as a coupling-constant integral over the 
density as a functional of the potential. Conditions are derived under which potential-ixmction&\ 
approximations are variational. Construction via this method and imposition of these conditions are 
shown to greatly improve the accuracy of the non-interacting kinetic energy needed for orbital-free 
Kohn-Sham calculations. 
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In the original form of density functional theory 
(DFT), suggested by Thomas [1] and Fermi [2] (TF) and 
made formally exact by Hohenberg and Kohn[3], the en- 
ergy of a many-body quantum system is minimized di- 
rectly as a functional of the density. Its modern in- 
carnation uses the Kohn-Sham (KS) scheme [4], which 
employs the orbitals of a fictitious non-interacting sys- 
tem. This brilliant idea means only a small fraction 
of the total energy need be approximated, and good 
approximations^, 6] have made DFT the popular tool 
it is today. DFT [7] is now ubiquitous in many scientific 
fields, including both materials and chemistry. 

Interest is rapidly reviving in finding an or&ita/- free [8] 
approach to DFT. The major bottleneck in modern calcu- 
lations is the solution of the KS equations, which can be 
avoided with a pure density functional for the kinetic en- 
ergy of non- interacting fcrmions, T s . The original TF ap- 
proximation is of exactly this type, but is far too inaccu- 
rate for modern applications. Despite decades of effort [9] , 
no generally applicable approximation for T s has been 
found, although material-dependent approximations [10] 
have been suggested, or approximations designed only 
for weakly-interacting systems [11]. 

However, Englert and Schwinger[12] pointed out that 
the potential is a more natural variable to use in deriv- 
ing approximations to quantum systems. In particular, 
semiclassical approximations begin with the classical mo- 
mentum, a local functional of the potential. TF theory 
is often derived first in terms of the potential, which is 
then eliminated in the final expressions, yielding an ex- 
plicit density functional. Exact potential functional the- 
ory (PFT) satisfies a variational principle with minimiza- 
tion over trial potentials [13, 14], yielding useful insight 
into the optimized effective potential method [13]. 

In the present work we go beyond those results by con- 
sidering explicit potential functional approximations to 
interacting and non- interacting systems of electrons; such 
approximations are presently being developed via a sys- 
tematic asymptotic expansion in terms of the potential, 



which has already been found in simple cases [15, 16]. The 
leading terms in a semiclassical expansion yield local ap- 
proximations to the energies, and the leading corrections 
greatly improve over the accuracy of local approxima- 
tions in a systematic and understandable way. Correc- 
tions to TF are relatively simple functionals of the po- 
tential, but far more subtle as functionals of the density. 
Such expansions are significantly more accurate and less 
problematic for the density itself rather than for the ki- 
netic energy density because the latter requires two spa- 
tial derivatives [16]. 

By minimizing over TV-particle wavefunctions that 
are antisymmetric, normalized, and have finite kinetic 
energy, we obtain the ground-state (gs) energy 
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as a functional of the potential, where T is the ki- 
netic energy operator, V ee the electron-electron repulsion, 
and V the one-body potential. We define the potential 
functional [13]: 



F[v] = <*„ | f + V ee | #„} 



(2) 



with tyy denoting the gs wavefunction of potential v(r). 
With F[v], the exact relation for the gs energy is 



E v = F[v] + / d 3 r n[v]{r) v{r) 



(3) 



in practice, requiring approximations to both F[v] and 
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F A [v] + d 3 r n A [v](r) v(r) , 



(4) 



where A denotes an approximation as a functional of the 
potential. We call this the direct approach. 

We show that (i) the universal functional, F[v], is de- 
termined entirely from knowledge of the density as a 
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functional of the potential, such that only one approxi- 
mation is required, namely n A [i;], (ii) the variational prin- 
ciple imposes a condition relating energy and density ap- 
proximations, (iii) a simple condition guarantees satisfac- 
tion of the variational principle, (iv) with an orbital-free 
approximation to the non-interacting density as a func- 
tional of the potential, the kinetic energy is automatically 
determined, i.e., there is no need for a separate approx- 
imation, and (v) satisfaction of the variational principle 
improves accuracy of approximations. 

We deduce an approximation to F from any n A [u](r) 
in the following way. Introduce a coupling constant in 
the one-body potential: 

t> A (r) = (l-A)« (r) + Af(r), (5) 

where vq(t) is some reference potential (possibly 0). In 
the context of TF theory this coupling was used in 
Ref. [17]. Then, using the Hellmann-Feynman theorem, 
-l 
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where Av(r) = v(r) — vq(f). Defining n[u](r) 
f Q dX n[v x ](r) and choosing vq(t) = 0, we obtain 

F[v] = / d 3 r{fi[v]{r) - n[v](r)} v{r) . 



(6) 



(7) 



This formula establishes that the universal functional is 
determined solely by the knowledge of the density as 
functional of the potential. Moreover, insertion of ?i A (r) 
on the right defines an associated approximate -F cc [n A [t/]], 
where cc denotes coupling constant. 

On the other hand, much of the accuracy of DFT cal- 
culations derives from the variational principle. In PFT, 
this yields 



£ A ' var = min (^F A [v} + J d 3 r n A [i]{r) v(r) \ , 



(8) 



with a possibly different value from Eq. (4) for a given 
pair of approximations. Experience suggests use of the 
variational principle improves results. The Euler equa- 
tion for the minimum is 

S 4ty = - frfV«(r')x A H(r',r), (9) 
6v(r) J 

where x A [v](r, r') = 5n A [v](r)/5v(r') denotes the density- 
density response function. If a pair of approximations 
{F A ,n A } satisfies Eq. (9) at v(r), then Eqs. (4) and (8) 
yield identical results, but this is not guaranteed a priori 
in approximate PFT. 

We next ask: Does a given i 7 ' cc [n A [v]] satisfy Eq. (9)? 
Taking the functional derivative of Eq. (7) yields Eq. (9) 
if, and only if, 
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(10) 



This condition is true in turn, if and only if, 

X >](r,r') = X >](r',r) . (n) 

The exact response function satisfies this relation, but 
an approximate functional n A [v] might not. This con- 
dition guarantees conservation of particle number under 
small changes in the potential, an elementary version of 
a conserving approximation[18]. 

A simple example illustrating these results is TF the- 
ory, considered as a potential functional. Then, 
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(12) 



where v s (r) = v(r) + v H (r), and the latter is the Hartree 
potential, determined self-consistently from 



W(r) = -47rn TF (r) 



(13) 



while [i is the chemical potential, determined by the nor- 
malization requirement that J d 3 r n(r) = N. Taking 
functional derivatives with fixed N, the usual TF energy 
expression [19] 

F TF [v] = T S TF M+ J d 3 rn TF [v](r) v{r)+U[n TF [v}] (14) 

satisfies the Euler condition when combined with 
Eq. (12), and F cc [n TF [v}} = F TF [v], where T S TF [v] = 
it 2 J d 3 r n TF [w](r) 3 /6 denotes the TF kinetic energy and 
U[n TF [v}] = Jd 3 r n TF [v](r) v H [n TF [v]](r)/2 the Hartree 
energy. 

In practice, the usefulness of these results for interact- 
ing electrons might be limited, as they require an ap- 
proximation to the interacting density as a functional of 
the one-body potential that is sufficiently accurate to be 
competitive with standard KS-DFT calculations, i.e., be- 
yond the accuracy of TF theory. Of much more practical 
use is their application to the non-interacting electrons 
of the KS scheme, which sit in the effective KS potential, 
which includes both a Hartree and (some approximate) 
exchange-correlation (XC) contribution: 



v s (r) = v(r) + 
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M(r') 



< c [<M](r). (15) 



For a given approximation to E xc , which determines v A c , 
this equation can be easily solved by standard iteration 
techniques, bypassing the need to solve the KS equations. 
A given n A [u s ] removes the need for solving any differen- 
tial equation in each iteration. 

However, once self-consistency is achieved, we need to 
extract the total energy of the interacting electronic sys- 
tem, for which we need the kinetic energy of the KS 
electrons. All our derivations apply equally to the non- 
interacting problem, so we deduce: 



T s [v] = / d 3 r {n s [v]{r) - n s [v]{r)} v(r) , 



(16) 



which is the analog of Eq. (7) for a system of non- 
interacting electrons in the external potential v(r) (which 
is called v s (r), when it is the KS potential of some inter- 
acting system) . This defines a kinetic energy approxima- 
tion determined solely by the density approximation: 

77KM] = / d 3 r{n$[v](r) -n>](r)} v(r) . (17) 



This is our main result for the non-interacting case. It 
eliminates the need for constructing approximations to 
the non- interacting kinetic energy T s . 

To illustrate the power of these results, we consider 
a simple example, a system of non-interacting, spinlcss 
fermions in a one-dimensional box. We choose vq(x) to 
be inside a box (0 < x < L), and oo outside. Then 
Av(x) = v(x) is some potential inside the box. For this 
case Eq. (17) (with a nonzero Vq) reduces to 



while the green curve is the approximation derived at 
great length in Ref. [16]. The small discontinuity at 
about x = 0.2 and 0.8 is where the approximation 
switches from a form that is asymptotically correct in 
the interior to one that is asymptotically correct near the 
walls. The error is shown in the bottom panel. Note that 
the approximation for t s (x) of Ref. [16] is already a con- 
siderable improvement over that used in Ref. [15]. The 
red and black curves result from Eq. (18). The black is 
exact, while the red uses the approximation for the den- 
sity in Ref. [15]. Their difference is plotted in the bottom 
panel, and is both locally and globally far smaller, and 
required no separate approximation for the kinetic energy 
density. 

The approximations of Refs. [15] and [16] were de- 
signed to be asymptotically exact as TV — >• oo, both for 
the density and the kinetic energy. In Table I we show 
errors compared to the exact result of from TF the- 
ory, the WKB approximation, Ref. [15], its improvement 



Ts>s Ml = E (N) + fdx {n*[v](x) ~ n£[v](x)} v{x) , in R- cf - [16], and T s cc in Eq. (18) with n£[v] of Ref. [15; 

Jo Even for TV = 1, T s cc is about two orders of magnitude 



(18) 

where E (N) = it 2 (TV 3 + 3/2 TV 2 + TV/2) /6 denotes the 
total non-interacting energy of TV spinless fermions in an 
infinite square well. 



TABLE I. Total non-interacting kinetic energy of N particles 
and its absolute error in TF, WKB, T* of Ref. [15] and [16], 
and T s cc of Eq. (18) for the potential of Fig. 1. 
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more accurate than WKB, and significantly more accu- 
rate than the direct approximations of Refs. [15] and [16]. 
As TV — » oo, T™ converges most rapidly. 
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FIG. 1. Exact (ex) and approximate (app) kinetic energy den- 
sities (above) of Eq. (18) (black, red) and of Eq. (19) (blue, 
green) with the approximation in Ref. [16], and their abso- 
lute errors (below) for one particle in v(x) = —5 sin 2 ("7r:r), 
< x < 1. (color online) 

In Fig. 1 we plot two distinct kinetic energy densi- 
ties, along with approximations to them, and the corre- 
sponding errors, for v(x) = — 5sin 2 (7ra) in a box of unit 
length. The blue curve is the exact kinetic energy density 
obtained from a traditional definition, 
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FIG. 2. Difference of total exact (ex) energy E and Eq (black) 
compared to the TF approximation (blue), T* of Ref. [16] 
(green) and T| c of Eq. (18) (red) for TV — 1 and trial potentials 
Av(x) = —D s'm 2 (irx). (color online) 
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We finally test the symmetry condition of Eq. (11). 
To do this, we perform a variational PFT calculation, 
implementing Eq. (18). We take a given external poten- 
tial (—5 sin 2 (7rx)), calculate exact gs wave functions with 
different potentials, and find their energy. Fig. 2 shows 
the results when the well depth D is varied. The exact 
result is a black curve, whose minimum occurs at D = 5. 
The blue curve is the result of TF theory, which satisfies 
the condition, but is not very accurate. The green curve 
is the approximation of Ref. [16], which, while more ac- 
curate, docs not minimize at the true potential. This 
demonstrates that the pair of approximations T^[v] and 
n* [v] given there do not yield the same answer variation- 
ally and directly. Note also that the direct evaluation 
(green curve at D = 5) is more accurate, when com- 
pared to the exact result (black curve at D = 5) than 
application of the variational principle (green curve at 
its mininum D ps 6.5). This is because these approxima- 
tions were derived semiclassically, and have uncontrolled 
errors. On the other hand, our coupling-constant approx- 
imation of Eq. (18) (red curve) is both far more accurate 
and minimizes at the true potential. 
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FIG. 3. Exact (ex) total energy (black) compared to the direct 
approximation of Ref. [16] (green) and T$ c of Eq. (18) (red) for 
N = 2 and trial potentials v(x) = — 5sin 2 (7ra;) — AZ?sin 2 (27rx) 
with the external potential at AD = 0. (color online) 

However, we further tested the coupling-constant ap- 
proximation by adding the perturbation — AD sin 2 (27ra;) 
with varying depth AD to the given external potential 
(— 5 sin 2 (7ra;)). For N = 2, as shown in Fig. 3, the min- 
imum is no longer at AD = 0, despite the high accu- 
racy. This demonstrates that Eq. (11) is not satisfied for 
all possible variations of the potential around the true 
one. However, the breakdown appears small and does not 
much diminish accuracy. This breakdown is due to the 
small normalization error in the semiclassical density [16]. 
If the error changes with the potential, the correspond- 
ing x A cannot be symmetric. We have also calculated 



the response function for that approximation and found 
non-symmetrical terms [19]. This shows the utility of our 
analysis: Accuracy is likely to be further improved if the 
result can be easily modified to enforce Eq. (11). 

An alternative to the coupling-constant method used 
here is the virial theorem, which yields the kinetic energy 
from the potential and density alone [20]. We recommend 
that version which has an origin-independent kinetic en- 
ergy density given in Ref. [21], satisfying 

V 2 ( s (r) = ^V{n(r)V« s (r)}, (20) 

where d denotes the dimension of space. While either the 
virial or the coupling-constant formulation can be applied 
to realistic systems, we use the coupling constant here be- 
cause our illustrations involve box-boundary conditions, 
which create complications for the virial theorem [22]. 

The coupling-constant formulation can be applied to 
realistic systems with potentials that vanish at large dis- 
tances, using t>o = for a reference. Then, to keep the 
particle number fixed, employ the device of putting the 
system in a large box whose size is taken to oo at the 
end of the calculation. Either expression has a great 
advantage over traditional density-functional approxima- 
tions, such as generalized gradient approximations. For 
such approximations, there is always an ambiguity in 
the energy densities; a term that integrates to zero over 
the entire space can always be added[23]. However our 
energies use an approximation to the density, which is 
uniquely determined for all r, and so can be used to iden- 
tify the relative contribution to the energy from different 
regions [19]. 

Which variation (coupling constant or virial) is most 
useful in practice awaits general-purpose approximations 
for the density as a functional of the potential for an ar- 
bitrary three-dimensional case. But at least it no longer 
awaits the corresponding kinetic energy approximations. 

We gratefully acknowledge funding from NSF under 
grant number CHE-0809859. 



[1] L. H. Thomas, Math. Proc. Camb. Phil. Soc. 23, 542 
(1927). 

[2] E. Fermi, Z. Phys. 48, 73 (1928). 

[3] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 
(1964). 

[4] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 
[5] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
[6] A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 
[7] C. Fiolhais, F. Nogueira, and M. Marques, editors, A 

Primer in Density Functional Theory, Springer- Verlag, 

NY, 2003. 

[8] Y. A. Wang and E. A. Carter, Orbital-free kinetic-energy 
density functional theory, in Theoretical Methods in Con- 



5 



densed Phase Chemistry, edited by S. D. Schwartz, chap- 
ter 5, page 117, Kluwer, Dordrecht, 2000. 
[9] R. M. Dreizler and E. K. U. Gross, Density Functional 
Theory: An Approach to the Quantum Many-Body Prob- 
lem, Springer- Verlag, 1990. 
[10] C. Huang and E. A. Carter, Phys. Rev. B 81, 045206 
(2010). 

[11] T. A. Wesolowski, Phys. Rev. A77, 012504 (2008). 
[12] B.-G. Englert and J. Schwinger, Phys. Rev. A 29, 2339 
(1984). 

[13] W. Yang, P. W. Ayers, and Q. Wu, Phys. Rev. Lett. 92, 
146404 (2004). 

[14] E. K. U. Gross and C. R. Proetto, J. Chem. Theory 

Comput. 5, 844 (2009). 
[15] P. Elliott, D. Lee, A. Cangi, and K. Burke, Phys. Rev. 



Lett. 100, 256406 (2008). 
[16] A. Cangi, D. Lee, P. Elliott, and K. Burke, Phys. Rev. 

B 81, 235128 (2010). 
[17] L. R. Pratt, G. G. Hoffman, and R. A. Harris, J. Chem. 

Phys. 88, 1818 (1988). 
[18] G. Baym and L. P. Kadanoff, Phys. Rev. 124, 287 (1961). 
[19] A. Cangi and K. Burke, long paper in prep (2011). 
[20] N. H. March and W. H. Young, Nucl. Phys. 12, 237 

(1959). 

[21] E. Sim, J. Larkin, K. Burke, and C. Bock, J. Chem. 

Phys. 118, 8140 (2003). 
[22] T. Cottrell and S. Paterson, Phil. Mag. 42, 391 (1951). 
[23] J. P. Perdew, V. N. Staroverov, J. Tao, and G. E. Scuse- 

ria, Phys. Rev. A 78, 052513 (2008). 



